Soliton "molecules": Robust clusters of light bullets 



Lucian-Cornel Crasovan*, Yaroslav V. Kartashov^, Dumitru Mihalache* , Lluis Torncr 
Institute of Photonic Sciences, Universitat Politecnica de Catalunya, ES 08034 Barcelona, Spain 



Yuri S. Kivshar 

Nonlinear Physics Group, Research School of Physical Sciences and Engineering, 
the Australian National University, Canberra, Australia 



m 
O 
o 

(N 
C 

a 
m 

(N 



Victor M. Perez-Garcfa 
Departamento de Matemdticas, E. T.S.I. Industriales, 
Universidad de Castilla-La Mancha, 13071 Ciudad Real, Spain 

We show how to generate robust self- sustained clusters of soliton bullets-spatiotemporal (optical 
or matter-wave) solitons. The clusters carry an orbital angular momentum being supported by 
competing nonlinearities. The "atoms" forming the "molecule" are fully three-dimensional solitons 
linked via a staircase-like macroscopic phase. Recent progress in generating atomic-molecular co- 
herent mixing in Bose-Einstein condensates might open potential scenarios for the experimental 
generation of these soliton molecules with matter- waves. 
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Solitons - non-spreading, self-sustained wave packets 
- are at the core of nonlinear science, thus they have 
been investigated and observed in a variety of settings 
during the last two decades 0. Today, one of the most 
challenging open frontiers of the field is the elucidation 
of complex soliton structures or "soliton molecules" to 
be constructed from a number of " atoms" , each being a 
fundamental soliton. However, multi-soliton structures 
found so far 0-0] tend to self-destroy through expan- 
sion or collapse, or at best exist as meta-stable states 
which break apart by small perturbations. Here we re- 
veal, for the first time to our knowledge, a physical mech- 
anism for generating clusters which are made of stable 
fully three-dimensional light bullets that propagate sta- 
bly over huge distances even in the presence of random 
perturbations in the initial conditions. The core of our 
approach is the use of two-color parametric solitons sup- 
ported by competing nonlinearities 0, > which allow 
both, to generate stable fully three-dimensional solitons 
and to reduce the soliton-soliton interactions and enhanc- 
ing the clusters robustness. The clusters are thus multi- 
colored, carry orbital angular momentum, and are linked 
via a staircase-like macroscopic phase distribution. Wc 
present the analysis for optical spatiotemporal solitons, 
but our findings are intended to stimulate further theo- 
retical and experimental research in the case of matter- 
waves in Bose-Einstein condensates Il3|-ll9|. 



Spatiotemporal optical solitons, the so-called "light 
bullets" (LBs), are self-sustained objects localized in all 
spatial dimensions and in time (f° r a recent 

overview see Ref.|2^|). They result from the simultane- 
ous balance of diffraction and dispersion by the medium 
nonlinearity, and a two-dimensional version has recently 
been generated in quadratic nonlinear media [2^ . On 
one hand, spatiotemporal solitons are challenging objects 



for fundamental research, as examples of stable localized 
objects in three-dimensional nonlinear fields are rare in 
physics. On the other hand, spatiotemporal solitons hold 
promise for potential applications in future ultrafast all- 
optical processing devices [2j]-|22j, where each soliton 
represents a bit of information and should be employed 
for digital operations. Mult i- channel all-optical soliton 
networks have been proposed based on the concept of 
soliton clusters [3(1 |. the structures carrying many inter- 
acting individual solitons, recently introduced for two- 
dimensional solitons in saturable nonlinear media [y] . 

Soliton clusters can be viewed as a nontrivial general- 
ization of " spinning" solitons (or doughnut-like vortices) 
[13- H3 and necklace-ring beams an d they also 

appear in the study of active nonlinear systems such as 
externally driven optical cavities j^j, 39]. But the soli- 
ton clusters investigated so far tend to be unstable or 
meta-stable under the action of small perturbations. We 
have recently shown in the case of two-dimensional spa- 
tial solitons, that the competition between quadratic and 
cubic nonlinearities reduces the strength of the soliton- 
soliton interactions, thus makin g sp atial soliton clusters 
more robust under propagation [4fJ. Here, we consider 
for the first time the case of clusters made of fully three- 
dimensional light bullets, and show that they propagate 
stably over huge distances even in the presence of random 
perturbations. 

We consider the propagation of two-color (fundamen- 
tal wave and second harmonic) LB "molecules" (see the 
sketch in Fig. 1) in a bulk dispersive medium with com- 
peting quadratic and cubic (Kerr) self-defocusing non- 
linearities. Under suitable conditions, the interaction be- 
tween a fundamental frequency (FF) signal and its second 
harmonic (SH), in the presence of the self-defocusing cu- 
bic nonlinearity, dispersion and diffraction in the (3+1)- 
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dimensional geometry, can be described by the reduced 
model Hl-Eil 
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Here, T, X, Y and Z are the normalized reduced time, 
transverse spatial coordinates, and propagation distance, 
U and v are envelopes of the FF and SH fields, a measures 
the strength of the defocusing cubic nonlincarity, and 
j3 is a phase mismatch between the FF and SH waves. 
Here * stands for the complex conjugate of a complex 
field. Equations JQ) assume different group-velocity dis- 
persion coefficients at the two frequencies, a being their 
ratio, and assumes that the temporal group- velocity mis- 
match between them has been compensated. Notice that 
Eqs.(l) correspond to the simplest model of light prop- 
agation in media with competing nonlinearities (e.g., it 
assumes a non-critical, type I, oo or ee wave interaction). 
In practice, the strength of each of the possible cross- 
phase-modulations depends critically on the crystalline 
symmetry of the particular material employed through 
the polarizations of the fields involved, hence the actual 
value of the relevant elements of the nonlinear suscepti- 
bility tensor. However, Eqs. (1) are expected to capture 
the essential physics behind the soliton cluster evolution. 
The interaction Hamiltonian of the system is: 
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is a conserved quantity during evolution. Its absolute 
and local minima correspond to stable and metastable 
configurations, respectively. 

Circular light-bullet necklaces were constructed as su- 
perposition of N fundamental spatiotemporal solitons 
with different phases such that the overall phase jump 
around the core is a multiple of 2tt (see Fig. 1). We thus 
have: 
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where uq , vq are the fundamental solitons at both fre- 
quencies, r n are the soliton locations, whereas the soli- 
ton phases at those points are <j> n = 2nirM/N and 2(j> n 



, respectively. Here M determines the full phase twist 
around the cluster and plays the role of a topological 
charge ("spin"). We have considered circular soliton ar- 
rays, i.e. equally spaced "atoms" displaced on a circle 
of radius Rq . First, by appropriate numerical tech- 
niques (a standard band-matrix algorithm to deal with 
the resulting two-point boundary- value problem) we have 
found the families of stationary solutions to Eqs. (1) - 
i.e. the fundamental (non-spinning) three-dimensional 
spatiotemporal solitons (uq ,vq). In fact, the stationary 
three-dimensional parametric soliton can be well approxi- 
mated by a super-Gaussian "ansatz" with suitable chosen 
amplitudes and widths for both the FF and SH fields. 

The parameters that play an important role in the dy- 
namics of the LB "molecules" are the necklace topolog- 
ical charge M, the number of "pearls" TV forming the 
cluster, the initial radius of the necklace Rq, the energy 
Elb of each constituent soliton, the wave- vector mis- 
match (3 and the strength of the defocusing cubic nonlin- 
earity a. In almost all of our calculations we have consid- 
ered the phase-matching of the interacting waves, taking 
thus (3 = 0. We have also set a = 1, assuming equal 
dispersions at both frequencies, and a = 0.2 as the dy- 
namical equations possess scaling properties with respect 
to a. By increasing the strength of the defocusing cubic 
nonlinearity one will slow down the interaction between 
the constituent " atoms" . Taking into account that the 
medium with competing nonlinearities supports stable 
spatiotemporal vortices (vortex tori) with unit topologi- 
cal charge when their energy exceeds a threshold |44j , we 
have studied in detail the dynamics of soliton " molecules" 
which have the total energy exceeding the corresponding 
stability threshold energy of the vortex soliton. Because 
the energy threshold for the existence of a stable vortex 
torus at a = 0.2 is E t h ~ 9120 , we have considered here 
clusters with N = 5 and = 6 solitons, each constituent 
having the energy Elb = 2100 , whereas for the cluster 
with N = 4 "atoms", the individual energy E LB = 2824 
was correspondingly higher. 

Firstly, we have studied the dependence of the cluster 
interaction Hamiltonian (or equivalently, the effective po- 
tential, defined as H(R )/H(oo j) on the initial radius Rq 
and on the necklace charge M. This quantity gives im- 
portant hints when looking for soliton bound states (see, 
e.g. Ref. 0, 0| for a detailed analysis). While the 
interaction Hamiltonian for the N — 4 clusters does not 
posses any minima whatever the topological charge is (see 
Fig. 2(a)), for N — 5 and N = 6, local minima of the 
Hamiltonian are present for charge M = 1. For N = 5 
the minimum is at i?o = 13.5, whereas for N = 6 the min- 
imum is at Rq — 12. In our simulations we have added 
normally distributed noise with zero mean and variance 
a noise = 0.1 to the input "molecules". Keeping M = 1, 
we have varied the initial cluster radius Rq around the 
minimum value given by the effective potential approach 
and have found a range of optimal values of the input 
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radii that minimize the mean radius oscillations of the 
soliton cluster. For N — 6 the value Rq = 12 lies in the 
optimal radius interval, whereas for N — 5, the value 
i?o = 16 assures small oscillations of the mean radius. 

In order to check the predictions given by the study 
of the effective potential, we have numerically solved 
Eqs. Q by using a finite-difference scheme based on 
a Cranck-Nicholson time discretization followed by a 
Newton-Picard iterative technique and the Gauss-Seidel 
method for solving the obtained system of equations. 
Transparent boundary conditions allowing the radiation 
to escape from the computation window have been imple- 
mented. We have monitorizcd the evolution of the mean 
radius of the cluster defined as: 

R(Z) = 1 J J J (X 2 +Y 2 +T 2 )^ 2 (\u\ 2 +\v\ 2 )dXdYdT 7 

(4) 

where E = JJJ {\u\ 2 + \v\ 2 )dXdYdT is the total en- 
ergy. If the initial radius i?o of the cluster is large, then 
the mean radius R(0) at the entrance of the nonlinear 
medium amounts to i?(0) « Rq. 

The evolution of clusters with N = 4 (Rq — 12) and 
N = 5 (Rq = 16) constituents is quite robust as shown 
in Fig. 3. The "molecules" undergo rotation and clean 
up the initial noise in the first stages of propagation. 
Our estimations for the angular velocity u> of the soliton 
clusters end up with ui = 0.0027 for the N = 4 cluster- 
shown in Fig. 3 and u) = 0.0014 for the N = 5 one. Thus, 
cluster rotations are observable after large propagation 
distances. Only after thousands of diffraction lengths 
a quasi-periodic shrinking and expansion followed by a 
decay into several unequal fragments is observed as seen 
in Fig. 4. Our soliton clusters are much more robust 
than the LB clusters in quadratic and cubic saturable 
materials that survive only a few diffraction lengths in 
the presence of initial random noise. 

The simulations with other necklace charges (M = 0, 
M = 2 and M = 3) for clusters composed of N = 5 
(R = 16) and N = 6 (R = 12), show that the LBs 
forming the M = "molecule" fuse in 100 - 150 propaga- 
tion units, whereas the soliton clusters with net charges 
M = 2 or M = 3 expand indefinitely. Detailed sim- 
ulations performed for the N = 6-light bullet clusters 
with M = 2 show that, by varying the initial clus- 
ter radius, the clusters formed with overlapping solitons 
(10 < Rq < 20) expand rapidly whereas the clusters 
built with well separated LBs (Rq > 22) have a mod- 
erate mean radius variation for a propagation distance 
over 600 diffraction lengths. Notice that for a typical 
diffraction length of a few mm, this corresponds to sev- 
eral meters, orders of magnitude larger than the feasible 
crystal lengths. Similar results were obtained for the non 
phase- matching case (/? ^ 0). 

We have also studied the influence of the initial phase 
distribution on the cluster dynamics by simulating the 



evolution of two configurations with identical intensity 
distributions but different phases. 

The first one, build as per Eq. 3, having a staircase- 
like phase, destroys finally, after thousands of diffrac- 
tion lengths, by splitting into two spatiotemporal solitons 
(Fig. 5(a)-(d)), while the second one, having a ramp- 
like phase mask (see the inset of Fig. 5(e)), develops 
into a vortex torus (Fig. 5(e)-(h)). Thus, we arrive at 
the conclusion that the key factor that impede the LB 
" molecule" with a staircase- like macroscopic phase to ex- 
cite a vortex soliton is the sequence of the phase edge- 
dislocations (see the inset in Fig. 5(a)) existing between 
the neighboring solitons which form the cluster. 

In summary, we have revealed a key physical mech- 
anism for creating truly three-dimensional light bullet 
clusters which survive under random perturbations of 
the initial conditions. We have generated such structures 
numerically for a nonlinear optical medium with compet- 
ing quadratic and cubic nonlinearities. The experimen- 
tal demonstration of the concept with light waves faces 
many important challenges, including the generation of 
single light bullet. This goal requires the elucidation of 
a material setting with high quadratic nonlinearity, suit- 
able group-velocity-dispersions and low one-photon and 
two-photon absorption at both FF and SH wavelengths, 
as well as small group-velocity-dispersion, together with 
adequate cubic nonlinearities. This is a formidable task, 
thus progress is being made slowly. In this context we 
would like to mention that it was shown recently that the 
strength of the cubic nonlinearity can be tuned by means 
of optical rectification ji^l even though at present the 
technique has been developed only for one-dimensional 
beams. 

However, although we showed the concept in the case 
of light waves, our study is important to other fields such 
as the physics of hybrid atomic- molecular Bose-Einstein 
condensates |4*^-|54j. Indeed, recent experiments demon- 
strated coherent mixing of atomic-molecular condensates 
[Hi) which under suitable conditions should be approxi- 
mately described by coupled equations for the macro- 
scopic wave functions similar to Eqs. Tak- 
ing into account that to date the experimental observa- 
tions of bright solitons in condensates are restricted to 
quasi-one dimensional geometries 0,0, the matter- 
wave analogue of our light bullet clusters would corre- 
spond to clusters of condensate drops existing without a 
trap. 
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FIGURE CAPTIONS 



Fig. 1 Cluster composed of six spatiotcmporal two-color 
solitons. The topological charge M of the soliton cluster 
is equal to one. (a) The fundamental frequency field and 
(b) the second harmonic field, (c) The phase distribution 
at fundamental frequency and (d) the phase distribution 
at the second harmonic. 

Fig. 2 Effective interaction potential versus initial clus- 
ter radius for (a) N = 4, (b) N = 5 and (c) N = 6 soliton 
clusters for different net topological charges. Typical os- 
cillations of the mean cluster radius of solitons clusters 
with "spin" M = 1 for (d) N = A, R = 12, E LB = 2824, 
(e) N = 5,R = 1 6, E LB = 2f00 and (f) N = 6, R Q = 12 
and E LB = 2100. 

Fig. 3 Stable evolution of soliton clusters with M = 1 
under superimposed input random noise. Shown are the 



contour plots for the N = 4 cluster: (a), Z = 0; (b), 
Z = 25; (c), Z = 50 and the contour plots for the N = 5 
cluster: (d), Z = 0; (e), Z = 25; (f), Z = 50. Only 
the (X, Y) slices at T — of the fundamental frequency 
component are shown; the second harmonic field exhibits 
a similar behavior. The other parameters are the same 
as in Fig. 2(d) for the N = 4 cluster and as in Fig. 2(c) 
for the N = 5 one. 

Fig. 4 Cluster evolution over long distances and the 
onset of symmetry breaking instability. Shown are the 
isosurfaces |it| = 1.1 for the N = 4 (a)-(d) and the N = 5 
cluster (e)-(h). The parameters are the same as in Figs. 
2(d) and 2(e). 

Fig. 5 Comparative evolution of two clusters with identi- 
cal intensity distributions but different phase masks. The 
net topological charge is the same (M = 1) in both situ- 
ations. Top panels: evolution of a six-soliton "molecule" 
with a step-like phase distribution; bottom panels: evo- 
lution of a six-soliton "molecule" with a ramp- like phase 
distribution. Shown are the isosurfaces |u| = 1.1. The 
insets in panels (a) and (e) show the initial phase mask. 
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